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Abstract. We present a model to explain the mass segregation and shallow 
mass functions observed in the central parts of dense and young starburst stel- 
lar clusters. The model assumes that the initial pre-stellar cores mass function 
resulting from the turbulent fragmentation of the proto-cluster cloud is signifi- 
cantly altered by the cores coalescence before they collapse to form stars. With 
appropriate, yet realistic parameters, this model based on the competition be- 
tween cores coalescence and collapse reproduces the mass spectra of the well 
studied Arches cluster. Namely, the slopes at the intermediate and high mass 
ends are reproduced, as well as the peculiar bump observed at 6 M©. This 
coalescence-collapse process occurs on short timescale of the order of one fourth 
the free fall time of the proto-cluster cloud (i.e., a few 10 4 years), suggesting that 
mass segregation in Arches and similar clusters is primordial. The best fitting 
model implies the total mass of the Arches cluster is 1.45 x 10 5 M©, which is 
slightly higher than the often quoted, but completeness affected, observational 
value of a few 10 4 M©. The derived star formation efficiency is ~ 30 percent 
which implies that the Arches cluster is likely to be gravitationally bound. 



1. Introduction 

Mass segregation is observed in many young stellar clusters, both compact and 
sparse, with the most massive stars being preferentially located in their central 
parts (e.g., Pandey et al. 1992; Malumuth & Heap 1994; Brandl et al. 1996; 
Hillenbrand & Hartmann 1998; Fisher et al. 1998; Figer et al. 1999b, Sagar 
et al. 2001; Stolte 2002; Le Duigou & Knodelseder 2002; Sirianni et al. 2002; 
Gouliermis et al. 2004; Sharma et al. 2007). Among young clusters, the class 
of very dense clusters, harboring large numbers of massive stars and known 
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as starburst clusters such as the Arches cluster, NGC 3603, the Quintuplet 
cluster, Westerlund 1 and 2 and R 136 offers a very challenging testbed for 
the theories of massive star formation not only because of the extreme mass 
segregation observed in their central parts, but also due to the fact that their 
radially integrated stellar mass function is observed to be top heavy (Moffat et 
al. 1994; Figer et al. 1999; Stolte et al. 2002; Stolte et al. 2005; Stolte et al. 
2006; Kim et al. 2006; Harayama et al. 2007; Brandner et al. this volume). 
Similar trends are observed for the super star cluster NGC 1705-1 (Sternberg 
1998). 

In particular, the structure of the young (age ~ 1 — 2 Myrs), and dense 
(~ 10 5 stars/pc 3 ) Arches cluster which is located at a projected distance of only 
25 pc from the Galactic center has been studied extensively over the past few 
years. Its mass function has been determined using a variety of space borne 
and ground based instruments in the near infrared such as the Hubble Space 
Telescope (HST) NICMOS camera (Figer et al. 1999a), the Hokupa AO system 
on Gemini (Stolte et al. 2002), the NAOS/CONICA camera on the VLT (Stolte 
et al. 2005), and the NIRC2 instrument on Keck (Kim et al. 2006). Star 
counting in the inner annulus of the cluster (i.e., between the center and 1 time 
the core radius which is ~ 0.2 pc) is affected by incompleteness effects due to 
severe crowding, and the outer parts might still be contaminated by field stars. 
Thus, an interesting and well resolved area to test theoretical models is the 
second annulus. The converging results of these observations yield a very flat 
mass function with a ~ 1 {dN/dM = M~ a ; The Salpeter value is a = 2.35, 
Salpeter 1955) in the central annulus of the cluster, a value of a ~ 1.7 — 1.9 in 
the second annulus (i.e., between 1 and 2 times the core radius), and a nearly 
Salpeter-like or slightly steeper exponent in the outer areas. 

The mass segregation and shallow IMFs observed in starburst clusters such 
as Arches have been interpreted as the result of dynamical processes (Kim et al. 
2006, Portegies-Zwart et al. 2007). Although dynamical processes will inevitably 
play a role in enhancing the fraction of massive stars in the central regions of 
a cluster, two points remain unclear: The first one is that the primordial IMF 
might be shallow in nature which leaves little room for the subsequent effects 
induced by dynamical friction. The second issue is that for young cluster as the 
Arches, it is not clear whether dynamical processes have time to play a significant 
role (see Freitag et al. this volume). Moreover, models based on mass segregation 
by dynamical processes do not reproduce, to date, some features of the Arches 
cluster mass function such as the peculiar bump at ~ 6 Mq. In this work, we 
discuss a model to explain the observed mass function of the Arches and similar 
starburst clusters based on the coalescence of pre-stellar cores (PSCs) in the 
proto-cluster cloud and their subsequent collapse into stars. 



2. The Coalescence-Collapse model 

The coalescence model is discussed in detail in Dib et al. (2007a). Here, its basic 
features are summarized. PSCs are embedded in an isothermal molecular cloud 
(MC) (temperature is T = 10 K), at different locations r from the cloud's center. 
Both the PSCs and the MC are axisymmetric (PSCs are initially spherical but 
are likely to quickly flatten as time evolves). The MC and the PSCs have density 
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Figure 1. Mass spectrum of Jeans unstable pre-stellar cores in shells of 
width=0.025 pc located at different distances from the cloud center (at 0, 1, 
2„ 5, and 10 R c ; R c = 0.2 pc is the proto-cluster cloud core radius), where j3 
(here (3 — 1.8) is the turbulent velocity power spectrum. Bottom: cumulative 
number of cores in the regions between [0 — R c ], [0, 2 R c ], and [R c — 2 R c ). 



profiles which are flat in their inner parts and follow r -2 and r -4 laws in their 
outer parts, respectively. However, due to their different locations in the MCs, 
PSCs of a given mass M will have smaller initial radii R p (r, M) in the inner 
parts of the cloud. As time advances, the radius of the PSC will decrease due 
to gravitational contraction. The PSC contracts on a timescale, t contjP , which is 
equal to a few times its free fall timescale, and can be parametrized as: 

/ 3 \ 1/2 

where v > 1 and p p is the radially averaged density of the PSC of mass M p , 
located at position r in the cloud. Thus, the time evolution of the radius of the 
PSC can be described by the following equation: 
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R p (r, M, t) = Rp(r, M) e -^/ tcont ^ . (2) 

Once the instantaneous radius of a PSC of mass M p , located at position 
r from the cloud's center is defined, it becomes possible to calculate its cross 
section for collision with PSCs of different masses. The cross section for the 
collision of a PSC of mass Mi and radius Ri with another of mass Mj and radius 
Rj and which accounts for the effect of gravitational focusing is given by: 



<j(Mi, Mj,r, t) = TT [Rp^r, M t ,t) + R pJ (r, Mj,t)] 

2G(Mj + Mj) 

2v 2 (R Pii (r,M h t)+ R~ (r, Mj,t)) 



1 + 



(3) 



As molecular clouds are unlikely to be virialized (e.g., Dib et al. 2007b), it is 
assumed that the cores relative velocity is not constant at a given position in the 
cloud and is rather better described by a Larson type relation v(r) = VQr(pc) a 
(Larson 1981; Vq = 1.1 km s _1 ), with a lower limit being the local thermal 
sound speed, which is uniform across the isothermal MC. As initial conditions 
for the PSCs mass distributions at different cloud radii, we adopt distributions 
that are the result of the gravo-turbulent fragmentation of the cloud, following 
the formulation given in Padoan & Nordlund (2002). In these models, the mass 
distribution of cores is given by the following function: 



N(r, M) dlog M = /„(r) M -3 ^ 4- ^ 



where 



P(Mj)dMj 



d log M, (4) 



P(Mj) dMj 



2 M 



jo M -3 



2^ 



j exp 



In Mj - A 



o-d 



dMj. 



(5) 



In Eq. [J fo(r) is the local normalization factor, and (3 is the exponent of 
the turbulent velocity power spectrum and is related to a by j3 = 2a + 1. The 
dependence on the local dynamical and thermodynamical conditions is hidden in 



the terms A and ad which are given by A = InMjQ—ln n' and ad = ln(l+J\A 7 
where Mjq is the Jeans mass at the mean density no, n' the number density in 
units of no, M the local Mach number, and 7 a numerical factor found through 
numerical simulations to be 7 = 0.5. 



3. Results 

With the initial conditions described in §. 2, we follow the time evolution of the 
PSCs mass spectrum by solving the following equation of N(r, M, t): 



dN(r,M,t) 

= 0.5 x rjir) 

dt 
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fAM 

/ N(r, m, t) N(r, M — m,t) a(m, M — m,r, t) v(r) dm 

N(r,m,t)a(m,M,r,t)v(r) dm, (6) 

where the first and second terms in the right hand side of Eq. [6] correspond to the 
rate of creation and destruction of a PSC of mass M, at location r, respectively 
(Nakano 1966; Shadmehri 2004). In Eq. El AM = M — M min , and ry(r) is a 
coefficient which represents the coalescence efficiency, with rj < 1. This efficiency 
can be the result of various physical processes which can affect the coalescence of 
PSCs. For simplicity, we shall assume that r\ is independent of position. In order 
to evaluate the transition from PSCs to stars, we compare, at each timestep, 
the local coalescence timescale to the local contraction timescale for PSCs of a 
given mass. The local coalescence timescale is t coa i(r, M) = l/w coa i(r, M) where 
Wcoal is the coalescence rate (Elmegreen & Shadmehri 2003): 



2 1 /2 u ( r ) 



mbin 



Wcoai(r, M) = " ' u y ]T (Ri + Rj) 2 

V s hell(r) 



2G(Mj + Mj] 
2v 2 {Ri + R 3j 



(7) 



where mbin is the number of mass bins, and V s heii is the volume of the shell of 
width dr located at distance r from the MC's center. The contraction timescale 
is given by Eq. [TJ Whenever the local contraction timescale is shorter than 
the local coalescence timescale, PSCs are collapsed into stars. When a PSC 
collapses to form a star, a fraction of its mass is re-injected into the proto- 
cluster cloud in the form of an outflow. This mass loss is accounted for in 
a purely phenomenological way by assuming that the mass of a star which is 
formed out of a PSC of mass M p is given by M*=?/> M p , where ip < 1. Matzner 
& McKee (2000) showed that ip can vary between 0.25 — 0.7 for stars in the mass 
range 0.5 — 2 Mq. As there is no evidence so far, for or against, whether this 
result holds at higher masses, it is assumed here that a similar fraction of the 
mass of a PSC is lost in the outflow independent of its mass. 

Fig. [2] displays the time evolution of the cumulative PSC populations in 
the region [R c q — 2 Rco]= [0.2-0.4] pc, which corresponds to the second annulus 
of the Arches cluster for a model with r\ = 0.5, v = 10, mass of the cloud 
M c = 5 x 10 5 Mq, the radius and core radius of the cloud R c = 5 and R c o = 0.2 
pc, respectively, an initial peak density of the PSCs p p0 = 10 7 cm 3 , a = 
0.37,V> = 0.58, and the fraction of the total cloud mass present in the cores 
e = 0.5. In the initial stages, the most massive PSCs, that have a larger cross 
section, coalesce faster than the less massive ones, essentially by capturing the 
numerous intermediate mass PSCs and causing a rapid flattening of the spectrum 
at the high mass end. By t ~ 0.07 t ffyC (i// )C = (Svr/^Gpc) 1 / 2 ~ 3 x 10 4 yr is the 
MC free fall timescale), a first generation of the smallest PSCs collapses to form 
stars. As time advances, more massive stars are formed in the shell (massive 
cores collapse later because of their lower average density) and in parallel the 
PSCs population decreases. By t ~ 0.1 iy/ iC the intermediate mass PSCs which 
constitutes the largest mass reservoir for coalescence collapses into stars (see 
Fig. [3]). At this time, the turnover in the PSCs mass spectrum is located at 
~ 8 — 10 Mq. Since the reservoir of intermediate mass objects is depleted, 
the remaining massive PSCs coalesce at a slower pace before they collapse. By 
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Figure 2. Time evolution of the pre-stellar cores mass function (left) and 
stellar mass function (right) in the region of the proto-cluster cloud between 
[Rc — 2 R c ] . The stellar mass function is compared to that of the Arches stellar 
cluster mass function (Kim. et al. 2006; open squares). Fits to the simulated 
IMF (bottom right figure) yield slopes of -2.04± 0.02 and -1.72± 0.01 in the 
mass range of [1-3] M and > 15 M , respectively, in very good agreement 
with the observations. Fits over-plotted to the data are shifted by one dex 
for the sake of clarity. 



t ~ 0.25 tffc, all PSCs of different masses in the shell have collapsed and formed 
stars. Because of mass loss, the stellar IMF is shifted to lower masses (bump 
shifted to ~ 6 Mq). In summary, the resulting IMF is not very different from the 
PSCs mass spectrum after the initial and rapid stage of strong coalescence until 
t ~ 0.01 tff c . In Fig. [21 over-plotted to the numerical result is the cumulative 
mass spectrum of the Arches cluster in the annulus of [0.2-0.4] pc (Kim et 
al. 2006). The coalescence-collapse model agrees better with the observations 
than models based on mass segregation by dynamical friction. In particular, 
the bump at ~ 6 M is reproduced. Fits to the stellar spectrum yield slopes 
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of a = -2.04 ± 0.02 and -1.72 ± 0.01 in the mass ranges of [1 - 3] M and 
> 15 M & , respectively, in very good agreement with observational values. With 
this set of parameters, the mass of the Arches clusters is found to be 1.45 x 10 5 
which is slightly higher than the often quoted, completeness affect, observational 
value of a few 10 4 M . This implies a star formation efficiency of (Mass of the 
cluster/Mass of the proto-cluster cloud)=(1.45 x 10 5 M / 5 x 10 5 M o )=0.29 
which again implies that the Arches is likely to be a gravitationally bound cluster 
according to the results of Geyer & Burkert (2001). 



1 F 




0.00 0.02 0.04 0.06 0.08 0.10 0.12 0.14 
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Figure 3. Time evolution of the ratio of the pre-stellar cores contraction 
timcscale (equal to few times the free fall time) to the coalescence timescale 
for cores of masses, 3, 10 and 12 M located at the center of the R c — 2 R c 
annulus in the proto-cluster cloud. Time is in units of the proto-cluster cloud 
free-fall time t// ;C ^3 x 10 4 years. 



4. Conclusions 

The origin of the shallow mass functions and mass segregation observed in young 
starburst stellar clusters is explained by means of a model based on the efficient 
coalescence of pre-stellar cores and their subsequent collapse to form stars. The 
coalescence of cores causes a strong flattening of the turbulence generated, pri- 
mordial core mass function. Once the reservoir of intermediate mass cores is 
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depleted, and the radii of cores are significantly reduced by gravitational con- 
traction, coalescence proceeds at a much slower pace until all cores in the whole 
mass range collapse to form stars. The entire process is fast, lasting for a frac- 
tion (~ 0.25) of the proto-cluster cloud free-fall time, i.e., a few 10 4 years. When 
applied to the Arches cluster, the model is able to reproduce the observed mass 
slopes of a ~ 2 and ~ 1.7 at the intermediate and high mass ranges, respec- 
tively, as well as the peculiar bump at ~ 6 M . Using this model, we estimate 
a total mass of Arches of 1.45 x 10 5 M , which is slightly higher than the of- 
ten quoted, but completeness affected, observational estimates of a few 10 4 M . 
The star formation efficiency in this model is ~ 30 percent which implies that 
the Arches cluster is gravitationally bound according to the results of Geyer &; 
Burkert (2001). 

Acknowledgments. S. D. would like to thank the conference organizers 
for the partial financial support which enabled him to attend the conference. 



References 

Brandl, B., Sams, B. J., Bertoldi., Eckart, A., Genzel, R., Drapatz, S., Hofman, R. 

Loewe, M., & Quirrenbach, A. 1996, ApJ, 466, 254 
Dib, S., Kim, J., & Shadmehri, M. 2007a, MNRAS, 381L, 40 

Dib. S., Kim, J., Vazquez-Semadeni, E., Burkert, A., & Shadmehri, M. 2007b, ApJ, 
661, 262 

Eisenhauer, F., Quirrenbach, A., Zinnecker, H., & Genzel, R. 1998, ApJ, 498, 278 
Elmegreen, B., & Shadmehri, M. 2003, MNRAS, 338, 817 

Figer, D. F., Kim, S. S., Morris, M., Serabyn, E., Rich, R. & McLean, I. S. 1999a, ApJ, 
525, 750 

Figer, D. F., McLean, I. S., & Morris, M. 1999b, ApJ, 514, 202 

Fisher, P., Pryor, C, Murray, S., Mateo, M., Richtler, T. 1998, AJ, 115, 592 

Geyer, M. & Burkert, A. 2001, MNRAS, 323, 988 

Harayama, Y., Eisenhauer, F., & Martins, F. 2007, ApJ, accepted, (|arXiv:0710.2882p 

Hillenbrand, L. A., & Hartmann, L. W. 1998, ApJ, 492, 540 

Kim, S. S., Figer, D. F., Kudritzki, R. P., & Najarro, F. 2006, ApJ, 653, L133 

Larson, R. B. 1981, MNRAS, 194, 809 

Le Duigou, J.-M. & Knodelseder, J. 2002, A&A, 392, 869 

Malumuth, E. M. & Heap, S. R. 1994, AJ, 107, 1054 

Matzner, C. D., & McKee, C. F. 2000, ApJ, 545, 364 

Moffat, A. F. J., Drissen, L. & Shara, M. M. 1994, ApJ, 436, 183 

Nakano, T. 1966, Prog. Theo. Phys., 36, 515 

Pandey, A. K., Mahra, H. S. & Sagar, R. 1992, Bull. Astr. Soc. India 1992, 20, 287 
Padoan, P. & Nordlund, A. 2002, ApJ, 576, 870 

Portegies-Zwart, S., Gaburov, E., Chen, H.-C, Atakan Gurkan, M. 2007, MNRAS, 
378, L29 

Sagar, R., Munari, U., & de Boer, K. S. 2001, 327, 23 
Salpeter, E. E. 1955, ApJ, 121, 161 
Shadmehri, M. 2004, MNRAS, 354, 373 

Sharma, S., Pandey, A. K., Ojha, D. K., Chen, W. P., Ghosh, S. K., Bhatt, B. C, 

Maheswar, G., & Sagar, R. 2007, MNRAS, 380, 1141 
Sirianni, M., Nota, A., De Marchi, G., Leitherer, C, & Clampin, M. 2002, ApJ, 579, 

275 

Stolte, A., Grebel, E. K., Brandner, W., & Figer, D. F. 2002 A&A, 394, 459 
Stolte, A., Brandner, W., Grebel, E. K., Lenzen, R., Lagrange, A.-M. 2005, ApJ, 628, 
L113 

Stolte, A., Brandner, W., Brandl, B., & Zinnecker, H. 2006, ApJ, 132, 253 



